"moths" <-
structure(list(meters = c(25, 37, 109, 10, 133, 26, 4, 3, 3, 
27, 16, 6, 17, 3, 5, 163, 10, 5, 13, 63, 4, 4, 33, 241, 18, 2, 
182, 48, 20, 3, 36, 233, 44, 35, 8, 55, 6, 90, 44, 21, 36), A = c(9, 
3, 7, 0, 9, 3, 0, 0, 0, 39, 7, 12, 6, 2, 1, 5, 2, 2, 23, 10, 
5, 6, 2, 4, 2, 3, 4, 3, 3, 1, 3, 6, 1, 9, 10, 2, 0, 6, 0, 0, 
9), P = c(8, 20, 9, 2, 1, 18, 5, 5, 2, 5, 16, 0, 0, 0, 0, 1, 
4, 0, 6, 12, 0, 5, 1, 1, 0, 1, 2, 3, 4, 0, 1, 3, 1, 0, 0, 0, 
2, 2, 4, 4, 1), habitat = structure(c(5, 7, 3, 3, 8, 2, 2, 2, 
2, 5, 7, 3, 3, 3, 3, 8, 4, 2, 5, 7, 3, 3, 3, 8, 4, 2, 8, 4, 6, 
6, 6, 8, 4, 6, 6, 6, 2, 4, 4, 1, 6), .Label = c("Bank", "Disturbed", 
"Lowerside", "NEsoak", "NWsoak", "SEsoak", "SWsoak", "Upperside"
), class = "factor")), .Names = c("meters", "A", "P", "habitat"
), class = "data.frame", row.names = c("1", "2", "3", "4", "5", 
"6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", 
"17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", 
"28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", 
"39", "40", "41"))
